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Quartessence is one of the alternatives to A-CDM that has lately attracted considerable interest. 
According to this unifying dark matter/energy scenario, the Universe evolved from an early non- 
relativistic matter-dominated phase to a more recent accelerated expansion phase, driven by a single 
fluid component. Recently, it has been shown that some problems of the quartessence model, such 
as the existence of instabilities and oscillations in the matter power spectrum, can be avoided if a 
specific type of intrinsic entropy perturbation is considered. In the present article we explore the role 
of skewness in constraining this non-adiabatic scenario. We show that non-adiabatic quartessence 
and quintessence have different signatures for the skewness of the density distribution on large scales 
' and suggest that this quantity might prove helpful to break possible degeneracies between them. 
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According to the current standard cosmological model, the dynamics of the universe would be dominated by two 
unknown components: dark-matter (DM), responsible for structure formation, and dark-energy (DE), that causes the 
accelerated expansion. Although there are several candidates for both DM and DE, there is still no evidence of either 
of them in laboratory physics. From the point of view of simplicity, it would be interesting to explore the possibility 
J> ' that a single component plays the role of both DE and DM, reducing from two to one the unknown constituents of 
OO \ the universe. A model that provides a single description of DE and DM through "unifying-dark-matter" or simply 
quartessence [l| has attracted a lot of interest recently. A prototype of this model is given by the quartessence 
J£j I Chaplygin model (QCM) @. 

Both the background and linear fluctuations were extensively studied for QCM, and were compared to observational 
data. The generalized Chaplygin gas (as quartessence) appears to be compatible with all available data regarding 
the expanskm history (see e.g. ref. |3| and refs. therein). For adiabatic perturbations, a linear analysis was done for 



the CMB jj] and LSS [5|. In this case, only QCM models close to the "ACDM limit" are allowed. Recently, it was 
shown that problems (pointed out in [Q), such as the existence of instabilities and oscillations in the matter power 
spectrum of QCM, can be avoided if a specific type of intrinsic entropy perturbation is considered. Such non-adiabatic 
O | model is consistent with the 2dF power spectrum for any value of the model parameters in the permitted interval, 
as long as the effective shape parameter assumes certain values Q. An "averaging problem" was also pointed out 
as a shortcoming of quartessence However, it is straightforward to show that the above mentioned non-adiabatic 
quartessence does not suffer this kind of problem . 

Thus, up to the present time, we can say that adiabatic quartessence is disfavored by the data, but it is not 
possible to distinguish non-adiabatic quartessence from concordance models like ACDM and quintessence using the 
previously considered observables. However, as we shall see, measurable differences in the predictions of these models 
appear clearly in the nonlinear regime, in particular in the skewness of the matter distribution in large scales. In 
our investigation we specifically consider three different quartessence models. All of them have the ACDM model 
as a limiting case for the background solution. The analysis of these three cases indicates that our result should be 
applicable to more generic quartessence models. 

While most studies of the nonlinear regime deal with the clumping of pressureless fluid (DM), in the case of 
quartessence it is imperative that one includes the effects of pressure. Accordingly, we apply, and somewhat extend 
to include relativistic pressure, a method for the computation of density cumulants developed in refs. [IfllTlllT2T|. 
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II. GRAVITATIONAL GROWTH FOR QUARTESSENCE IN THE SPHERICAL COLLAPSE 

APPROXIMATION 

In our approach we consider the following newtonian-like equations fl3l IbH ] 

V^ = 4ttG(p + 3P), (I) 

— j +{u-V r )u = -V r( j)-{ P + P)- 1 V r P , (2) 
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V r • (pu) + PV r -u = 0. (3) 



Equations @ and © are, respectively, the Poisson, Euler and energy conservation equations, where relativistic 
effects of pressure (inertia and active gravitational mass) have been included. In these equations p, P, u and <fi stand 
for, the energy density, pressure, velocity field, and gravitational potential of the cosmic fluid, and we take c = 1. In the 
linear regime, for c e ff = (see bellow), system (11131 gives exactly the same equations as the relativistic perturbation 
theory in a particular gauge |l4j . Also the (nonlinear) energy conservation and the Raychadhuri equations derived 
from this system are formally identical to the general relativistic ones |l5j in this case, which provides a motivation 
for the above system. 

It is useful to split the dynamical variables into their background and inhomogeneus parts, i.e. we write: p = 
p + 5 p = p(l + 5) , P = P + SP, (p = (f> + ip, and u = Hr + v. Here H — a" 1 da/ dt is the Hubble parameter, a stands 
for the scale factor, and the overbar "~" denotes background (average) quantities. Introducing comoving coordinates 
x = f/a, neglecting shear and vorticity, and taking into account the background equations (assuming critical density), 
we obtain, after some algebra, the following differential equation for the density contrast 5, 
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In the above equation the prime symbol denotes differentiation with respect to r/ = In (a), H. = aH, w = P/p and 
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In the linear approximation, in the special case in which perturbations arc adiabatic we have 



c 2 = P'/p 1 — w — w' /3(l + w). However, as discussed in 0, in this case, the right hand side of 10} gives rise to 
oscillations and instabilities in the mass power spectrum that render the model unacceptable. This problem can be 
circumvented if non-adiabatic perturbations such that c? e ^ = 0, are considered 0- In the following we assume this is 
the case such that Eq. Q simplifies to, 
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To study the weakly nonlinear regime of structure formation and compute the higher order moments of the density 
distribution, it is useful to expand d as EH 113 : 



— T d o - 
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where do is a small perturbation. Using the above expansion we obtain for the linear (Di) and second order (-D2) 
factors the following differential equations, 
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FIG. 1: Contours of constant skewness in the (M 4 /po,a) plane for QCM (solid curves). We also plot curves of constant Q,m* 
(dashed curves) and constant present value of the equation of state wo (doted curves). 

Analogously, higher order modes can be obtained recursively by using the solutions of the differential equations for 
the lower order terms. Of special interest is the second-order equation. If we start with Gaussian initial conditions, 
it is associated with the emergence of non-Gaussian features in the matter density field. Further, D 2 can be related 
to the skewness of the cosmic field miii. In this case, the unsmoothed skewness is given by S3 = 3D2/ D\. 
We now consider the following quartessence models: 

P 

P 

P 

When a = 0, all the models above have ACDM as a limiting case for the background. At first and higher order in 
perturbation theory, however, these quartessence models with a = 0, have distinct behavior as compared to that of 
dark-matter in ACDM. 

In Figs. 1, 2 and 3 we show contours of constant skewness values at present time for the above models. In our 
numerical computation, we consider that at z = 10 3 , the growing modes D\ and D2 assume the Einstein-de Sitter 
behavior: D\ oc a, D2 oc a 2 , and such that S3 = 34/7 at that time. We also include in the figures two curves 
with constant effective matter density parameter (f2^ = 0.2 and 0.4). Here, = lim a ^o p/poa 3 , where po is 

the present value of the background energy density. We should expect the region between these two curves to be, 
roughly, the one allowed by current observational data. For instance, in Q it was shown that, for QCM, constraints 
from cluster x-ray data (that are the most restrictive ones), essentially correspond to these curves. The same holds 
for the other quartessence models. From the figures it is clear that, for all the considered models, the skewness in 
these regions assumes values between S3 r* 13 — 20. This strongly contrasts with what is expected in ACDM and 
quintessence, where one obtains S3 « 5, weakly sensitive to the cosmological parameters [ill 120L l2ll l22l| . 

At first sight this difference could be interpreted as an indication that quartessence models are inconsistent with 
large-scale skewness measurements [24|. However, care should be taken when analyzing this issue. Whereas in the 
discussion above baryons were neglected, measurements of skewness from large-scale galaxy distribution, are based on 
counting luminous objects, not the dark component. Eqs. (l)-(5) can be easily generalized to include baryons [23l |. 
The main outcome is that the quartessence skewness is not substantially affected by the presence of a small amount 
(flbo ~ 0.04) of baryons, but on the other hand, the baryonic skewness (S3b) is nearly constant with redshift, i.e. 
Ssb — 34/7. Therefore, again baryons behave differently from the dark component as in the adiabatic quartessence 



(Chaplygin Quartessence), (9) 
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(Logarithmic Quartessence), (10) 



-M 4 exp ( —j- ) (Exponential Quartessence). (11) 
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FIG. 2: Contours of constant skewness in the (M 4 /pn, a) plane for logarhitmic quartessence (solid curves). We also plot curves 
of constant fi^ (dashed curves) and constant present value of the equation of state wo (doted curves). 

power-spectrum case ||. We remark that this holds even for a = 0. 

Further investigation is necessary to clarify to what extent do skewness measurements from galaxy distribution 
constrain quartessence models. Although more challenging to observe, a potentially powerful probe is the lensing 
(convergence) skewness |25|, that is sensitive to both baryons and quartessence. According to [2(| current lensing 
observations are still too noisy to allow strong constraints on cosmological parameters. On the other hand, as discussed 
above, quartessence predictions for the skewness are quite different from ACDM models and its variants. Thus, we 
expect that present and upcoming data might discriminate among these two classes of models. 

In the present work, we have not quantitatively compared quartessence predictions with observations. Our goal here 
is more modest; we essentially use Figs. 1, 2 and 3 to stress out the important fact that non-adiabatic quartessence and 
concordance models like ACDM and quintessence could be observationally distinguished by skewness measurements. 
We leave this investigation for a future work [2i| . Finally, we remark that our results do not apply straightforwardly 
to quartessence models with a scalar field with, for instance, a non-canonical kinetic term p7^. but these models have 
yet to be tested against background and power spectrum data. 
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